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In order to get an estimate of the homogeneity of the distribution of matter in a fast hadron, we 
compute the correlation of the saturation scales between different impact parameters. We find that 
these correlations are quite strong: The saturation scale is nearly uniform in a wide domain around 
each point in impact-parameter space. We provide analytical expressions for the correlations, 
which are supported by numerical simulations. Although the numerical calculations are done for 
specific saturation models which are obtained from QCD after drastic simplifications, we expect 
our analytical formulas to be correct for full QCD in asymptotic limits, since their derivation 
requires only a few general assumptions. 
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1. Introduction 

Phenomenological models for the very high-energy regime of QCD where saturation effects 
become important (see Ref. [|l|] for a recent review) are usually built on some parametrization of 
the elastic dipole-hadron scattering amplitude T{y,r,b) which is a function of the rapidity y of 
the scattering,' of the size r of the dipole, and of the impact parameter b. This amplitude is then 
related to the observables through appropriate convolutions with the wave functions that describe 
the incoming objects. In the simplest of these models, due to Golec-Biemat and Wiisthoff [^1, the 
dipole amplitude is assumed to have the form 

T{r,y,b) = l-£.-'-'e.?(>''^)/4, (1.1) 

where the momentum scale Qs, called the saturation momentum, is parametrized as 

Q^,{y,b) = 1 GeV^ x 0(7? - b)e^^y-y°l (1.2) 

The constants R, X and 3^0 are determined from a fit to the inclusive deep-inelastic scattering data. 
The spatial distribution of matter in the plane transverse to the collision axis is encoded in the b- 
dependence of the saturation momentum. The 0-function used by Golec-Biernat and Wiisthoff is 
sometimes changed to a smoother distribution in such a way that the model be also able to describe 
semi-inclusive diffractive data. In any case, the fluctuations between different points in transverse 
space are completely neglected in all these models. Note that this may not be a problem for standard 
phenomenology since most of the observables in deep-inelastic scattering probe one single point 
in impact-parameter space in each event. But clearly, independently of phenomenology, we would 
like to understand better how the matter is distributed in a fast hadron. 

We shall first explain why fluctuations of the parton densities are expected between different 
impact parameters, then we shall provide a heuristical discussion of the form of these fluctuations, 
for which we have been able to write a parameter-free formula valid in some asymptotic limit. 

2. Picture of a fast hadron/nucleus 

Let us consider a fast hadron or nucleus probed by a color dipole of size r (which may be 
seen as a component of a virtual photon of virtuality 2 ~ 1 /r) at very high rapidity y. We go to 
a frame in which the probing dipole is almost at rest and we require that the impact parameter be 
some fixed b (see Fig. |l|). The scattering probability T is roughly proportional to the local density 
n of partons in the corresponding phase-space cell: T{r,y,b) ~ ajn{r,y,b). It proves useful to 
see r as a probability of interaction between the dipole and a fixed configuration of partons: The 
physical amplitude measured in experiments is then T averaged over events (i.e. over the partonic 
configurations; see Ref. [||] for a review). If the rapidity is high enough, we know that at each b, 
T has the shape of a front connecting 1 (black or saturated regime) for r ^ 1 / [y, b) to (color 
transparent or dilute regime) for r <C 1/ Qs{y,b). The saturation momentum Qs{y,b) determines the 
transition. It grows exponentially with y, which means that the position of the wave front moves 
linearly along the axis log(l/r^) when the rapidity increases. It was first conjectured ||3|] and then 

'Throughout our discussion, y is actually the rapidity multiplied by the factor a — asNc/n. 
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Figure 1: Schematic picture of a fast-moving hadron (left; the colored disks represent partons) and scattering 
amplitude as a function of the size r of the probing dipole at two impact parameters ^71,^2 for a fixed rapidity 
y (right). 



checked numerically [Q] that to a good approximation, at each point in impact-parameter space, the 
only variable relevant to the evolution of the amplitude T with the rapidity y is log(l/r^). 

In the regions in which T '-^ aj, few partons are probed, hence the further evolution of T with 
the rapidity y is stochastic. If T aj instead, many partons populate that phase-space region, and 
the evolution of T is of deterministic nature: A mean-field approximation of the QCD evolution 
can be taken. Thus for values of r of the order of the inverse saturation momentum, T has the 
shape of a smooth (deterministic) curve traveling towards smaller values of r. However, because 
of the fluctuations in the tail of the front, Qg is a stochastic variable for the rapidity evolution. 
Fluctuations in the dilute region of phase space propagate towards the dense region and affect the 
saturation momentum typically after an additional evolution over Ay ~ log^(l/ aj) units of rapidity. 
They result in a random diffusion of = log of variance {pj)c ~ Dy, where D can be computed 
from QCD 

These fluctuations determine a dispersion of from event to event. But stochasticity is also 
expected to manifest itself by differentiating the points say bi and b2 in impact-parameter space, 
creating a dispersion of Ps in the transverse plane. In order to characterize these fluctuations, we 
shall now compute the correlator = {{Ps{bi) — Pi (^2))^) at fixed y. 



3. How correlations may occur: heuristic discussion and analytical formulas 

Let us examine how correlations between two points in transverse space bi and b2 may build 
up. We define Ab = \b2 — b\\. If M) < \/Qs (Qs is the saturation momentum at either bi or b2), 
then obviously Qs{bi) = Qs{b2) and af2 = 0. If Ab > l/Qs instead, then the evolution around the 
impact parameter bi can influence the evolution around b2 only if a parton at bi splits into a parton 
of size of the order of Ab. But the saturation of the density of partons of sizes larger than l/Qs 
disfavors such splittings. Hence we may think that the evolutions decouple as soon as the saturation 
radius l/Qs becomes smaller than Ab. Assume that this happens at rapidity jo- Then for y < yo, 
0^2 = 0, and for y > yo, 0^2 — {Pi)c + (pI)^ ~ 20(3^ — ^0)- One may fix the rapidity y and vary 
the distance Ab instead: Then 0^2 ~ '^D\og{AbQs) / x' {Jq) for \og{AbQs) > (see the dotted hne 
in Fig. |2|), which suggests that the characteristic distance scale for the correlations in the transverse 
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Figure 2: Sketch of the correlations as a function of the logarithm of the distance Ab= \bi — in impact- 
parameter space scaled by l/Qs- For AbQ^ > 1, the points bi and b2 are statistically independent. The 
dotted line represents what one would naively expect if fluctuations affected the saturation scale as soon as 
^bQs > 1. (£> is the diffusion coefficient of p.y for a single front, namely D ~ {p]')c/y)- The continuous 
line takes into account the delay induced by the propagation of the fluctuations, which results in an effective 
persistence of the correlations. 



plane is l/Qs- x{7o) is a particular eigenvalue of the BFKL kernel and x'iYo) the asymptotic 
rate of change of (p^ ) with the rapidity [^]. 

However, this is not yet the correct answer. Indeed, as recalled before, for fluctuations to 
be able to differentiate bi and b2, ~ log^(l/a,^) extra units of rapidity are needed after the 
rapidity yo at which AbQs{yo) = 1- Hence the effective decoupling of the saturation momenta is 
expected later in rapidity, or for larger distances Ab. The correlations would persist over distances 
Ab ~ f.fiog'(V«r) /g^ (see the full hne in the sketch of Fig. |). 

Extending the phenomenological theory for stochastic fronts developed in Ref. |^], we are able 
to fully compute the correlator 0^2 = {{Ps{b) — Ps{b + Ab)Y). One way of writing the result is 



2 _ 2ti^ 

~ 3yglog(l/«2)yexp|-'''^^'''''°'N'i'^'^°;'°^'^''"-' 



dq 



(3.1) 



where 1^4 is a particular Jacobi theta function. The interesting limiting behaviors read 



3Z'(0 0) log' 
4 



X"{ro)iog{^hQ,) 



exp 



;t:'(yo)log^(lM^) 

2r^z"(ro)iog(Afce,: 



forlog(A/7a)»log'(l/a.') 

for log(l/a2) < log(A/7e,) < log2(l/a2). 

(3.2) 



Comparing the expression of af2 in the large Ab limit to the variance Dy of p^, we find that 0^2 is 
actually equal to 2Dlog{AbQs) / x' (Yo) for large log{AbQs). From the second limiting expression, 
it is obvious that (7^2 is close to zero forlog(AZ7Gv)<log^(l/«,2)- 

In order to check these expressions, we performed numerical simulations of models which 
possess the main characteristics of the QCD evolution while being simple enough to allow for 
robust Monte Carlo simulations (see Ref. [^] for details). We found a perfect matching with the 
parameter-free analytical result (3.1) in the limit log(l/a^^) ^ 1. For larger and more realistic 
values of a^ , the persistence of the correlations is still seen in the numerical simulations, but some 
parameters should be modified in the analytical expressions and tuned to account for our lack of 
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Figure 3: Comparison of a numerical Monte Carlo simulation and our analytical formula. The constant in 
the parameter L (see Ref. for the definitions of L and 5), which should be equal to log(l / aj) / 70 for very 
small Us, has been shifted by a phenomenological constant. Once this is done, we get a very good agreement 
between the two calculations. 

understanding of subleading corrections important for finite log(l/ aj). We show such a calculation 
for = 0.1 in Fig. ||, compared to a variant of Eq. (^). 

4. Conclusion and outlook 

The main result of our work is that the characteristic distance of the correlations in the trans- 
verse plane is not l/Qs as one would naively expect, but rather exp [clog^(l/a^^)] /Qs (c being a 
known constant), which is parametrically much larger than l/Qs. Our results are valid for large 
log(l /ccj), and for distances Ab much smaller than the typical confinement scale 1 /Aqq^. 

The goal of our work was to understand the fundamentals of the QCD dynamics in transverse 
space, without thinking a priori of any application to phenomenology. Let us however note that 
recently, a diffractive deep-inelastic scattering observable was proposed that would directly probe 
the correlations which we have computed [0]. (A calculation of these correlations in the framework 
of the B-JIMWLK formalism ^ which a priori neglects the fluctuations discussed in this paper 
was also performed.) Also, these correlations may play an important role in heavy-ion collisions. 
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